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In this work we study the Casimir effect with three-dimensional topological insulators including 
the effects of temperature and uniaxial anisotropy. Although precise experimental values for the 
optical properties of these materials are yet to be established, qualitative analysis is still possible. We 
find qualitatively that the reported repulsive behavior and the equilibrium point are robust features 
of the system, and are favored by low temperatures and the enhancement of the optical response 
parallel to the optical axis. The dependence of the equilibrium point with temperature and with the 
topological magnetoelectric polarizability characteristic of three-dimensional topological insulators 
is also discussed. 
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I. INTRODUCTION 

Since H.B.G. Casimir wrote his seminal paper— propos- 
ing that two uncharged metallic plates should attract in 
vacuum due to quantum fluctuations still little is known 
of how to revert the sign of the force and obtain repulsion 
or levitation. The generalization of the theory to include 
dielectric bodies^ has made it possible to achieve, very 
recently, repulsion with test bodies immersed in a di- 
electric fluid^. However, achieving repulsion in vacuum, 
which would have important consequences for device 
applications^'^, is still to be experimentally tested. The- 
oretically, there have been several proposals to achieve 
such situation. The first proposal considered magnetic 
versus non-magnetic situations^ and more recent pro- 
posals include the use of metamateriala^r— and geome- 
try to induce repulsive scenarios!*^. Importantly, it has 
also been shown that for general dielectric-dielectric sit- 
uations, there is no repulsive behavior in vacuum^^. In a 
recent paper—, two of us suggested that repulsive behav- 
ior could be achieved by using topological insulators^"— . 
The topological nature of these materialsA2ii§ provides a 
magnetoelectric term modifying the response to an ex- 
ternal electromagnetic field, which can then cause the 
Casimir effect to reverse sign at short distances and re- 
main attractive at large distances^^. Consequently, at an 
intermediate distance d^q , the Casimir force vanishes and 
a stable configuration is possible. 

In this work, we aim to include the effect of temperature 
and uniaxial anisotropy to this behavior, two relevant fac- 
tors in experimental situations. It is well known by know 
that anisotropy can modify the Casimir force by chang- 
ing the form of the refiection coefficients^, and has been 
studied in the context of repulsive interactions^, even in- 
cluding the effect of temperature^^. In the present paper 
we address the question of whether uniaxial anisotropy 
and temperature can enhance or destroy the repulsive 
behavior in the context of topological insulators, and in 
what cases can this occur. To achieve this goal we de- 



rive the refiection coefficients for a generic anisotropic 
topological insulator, a result which to our knowledge is 
absent in the topological insulator literature. With the 
help of these, we find that low temperatures favor the 
repulsive behavior and we determine that enhancing uni- 
axial anisotropy in the direction parallel to the optical 
axis also acts to increase repulsion. 

To our knowledge, precise optical characterization of 
these materials is still lacking, and thus it remains an 
open question if real materials are to show repulsion at 
these distances. Nevertheless, throughout the paper, we 
use choices of parameters which can be relevant for topo- 
logical insulators such as TlBiSe2^^^ (see also concluding 
paragraph in Ref.^^) which are mainly low dielectric per- 
meability at zero frequency e(0) and small magnetoelec- 
tric polarizability 6. We do not restrict the study only to 
these values but explore other regions of parameters for 
completeness. 

The work is structured as follows. Section [ll] is devoted 
to review the formalism and notation on which we rely 
on. To do so, we will briefiy review the isotropic case at 
zero temperature (T = 0) discussed in Ref J^, to which 
we will add the dependence of the minimum with the 
topological magnetoelectric polarizability, not discussed 
in Ref . . Details of the proof of the existence of repulsion 
are left for appendix [Xj In section IIIII we will study the 
effect of temperature on the Casimir energy as a function 
of distance, both analytically and numerically. We will 
also discuss the position of the minimum as a function 
of temperature and topological magnetoelectric polariz- 
ability 9. Details of the analytical results can be found 
in appendix [BJ In section IIVI we will discuss the effect 
of uniaxial anisotropy on the Casimir energy at T = 0. 
First, we will explicitly write the Fresnel coefficients for 
this case, leaving the details of the derivation for ap- 
pendix [Cj In the remaining part of this section we will 
consider different values for the relevant parameters and 
discuss their effect on the position of the minimum. Fi- 
nally, in section |V] we will present the results for the case 
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of uniaxially anisotropic topological insulator at T ^ to 
end with section IVTl which includes a summary of results 
and some concluding remarks regarding real systems. 
We emphasize that the numerical results presented in this 
paper explore regions of parameters without any restrict- 
ing assumption other than positive energy condition in- 
herent to magnetoelectric couplings, discussed in section 
IIIII B. On the other hand, in order to derive analytical 
results, some feasible assumptions are to be considered, 
such as small magnetoelectric coupling, all of which ap- 
ply only to the results in appendix [Bl and hence will be 
discussed in the mentioned appendix. 




II. CASIMIR-LIFSHITZ FORCE BETWEEN 
TOPOLOGICAL INSULATORS AT T = 

It is well known from different approaches, that the 
Casimir energy density stored by two dielectric parallel 
plates (labeled 1 and 2) of area A, separated by a distance 
d, at zero temperature (T = 0) can be written as^'^ i^^'^'^ : 



Ah 

(1) 

where fcs = \J^f\ + C^/^^ the wave vector perpendicu- 
lar to the plates, fc|| is the vector parallel to the plates 
and ^ is the imaginary frequency defined as lo — i^. The 
matrices Ri,2 are 2x2 reflection matrices containing the 
Fresnel coefficients such that: 
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i?,,,(ie,k||) i?.,p(ie,k||) 
■Rp,s(«^,k||) i?p,p(i5,k||) 
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where Ri j describes the reflection amplitude of an in- 
cident wave with polarization i which is reflected with 
polarization j. The label s (p), equivalent to TE (TM) 
modes, describes parallel (perpendicular) polarization of 
the electric field with respect to the plane of incidence. 
The Casimir force per unit area on the plates is obtained 
by differentiating expression ([T]) : F = —ddEc{d). A nega- 
tive (positive) force, or equivalently a positive (negative) 
slope of Ec{d), corresponds to attraction (repulsion) of 
the plates. 

In order to calculate the Casimir energy one needs to 
compute the reflection coefficients. These coefficients re- 
late the amplitude of the incident and refiected electric 
fields for different polarizations. Since the normal com- 
ponents of D and B and tangential components of E and 
H must be continuous along the interface, one can obtain 
the reflection amplitudes by solving Maxwell's equations 
at each side and then imposing the mentioned matching 
conditions for the fields (see below). 



A. Topological insulators: The model 

Following RefJ^ we will consider a Casimir system 
composed by two topological insulator plates separated 



FIG. 1: (color online)Topological insulating plates separated 
by a distance d covered with a thin ferromagnetic layer of 
thickness I « d (not to scale). Changing the sign of the 
magnetization of one of the plates (parallel or antiparallel 
to the surface normal) is equivalent to changing the sign of 
(either) 9i^2- 



by a distance d, shown in FiglT] The electromagnetic re- 
sponse of three dimensional topological insulators, which 
determines the refiection coefficients, is governed by the 
lagrangean: 



C^Co + Cg^{E-D + B-H) 



(3) 



with the difference that now the constitutive relations for 
the electric displacement field D and the magnetic field 
H are given by: 



Di - e,jEj 
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H, = {^,)TlB,-—E,, 



(4) 
(5) 



The first term of each constitutive equation is the 
usual electromagnetic term defined in terms of the di- 
electric function t{uj) and magnetic susceptibility /i(w), 
giving rise to the ordinary electromagnetic term in the 
lagrangean ^ which we label Cq. We restrict ourselves 
to non magnetic topological insulators so = 1. The 
second term in each constitutive equation is a non trivial 
axionio^i or topological magnetoelectric terroiiii^ which 
gives rise to a topological Ce term in the lagrangean. 
Time reversal symmetry indicates that = 0, tt (mod 
2n) being — tt the case for topological insulators and 
= Q the case for trivial insulators. 

When the boundary is included, the action corresponding 
to the lagrangean ([3]) , S'o -f- 5*9 , is a fair description of the 
topological insulator only when a time reversal breaking 
perturbation is induced on the surface to gap the surface 
statesii. Along the lines of similar situations^iiiiS^, in 
this work we consider that the time reversal perturbation 
is a magnetic coating of small thickness I « d {d will be 
of the order of /im) which gaps the surface states. In the 
described situation, 9 is quantized in odd integer values 
of TT such that: 
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e = {2n+l)-K, (6) 

where n S Z, determined by the nature of the coating, 
but independent of the absolute value of the magneti- 
zation of the coating. Positive or negative values of 6 
are related to different signs of the magnetization on the 
surface? - , which we consider is perpendicular to the plane 
of the plates. Each topological insulator plate is charac- 
terized by the frequency dependent dielectric function 
e{io) and the quantized magnetoelectric term 6. Being 
a topological contribution, 9 is defined in the bulk as a 
constant whenever the bulk Brillouin zone is definecU^. 
Nevertheless, at the surface, the effective action must be 
valid as long as the condition / << d is fulfilled, where 
magnetic effects from the covering magnetic layer could 
renormalize its value, and the system should be treated 
as a layered system. We henceforth assume that this is 
not the case and consider only the region where I « d 
(nm << /im in practical terms). 

The parasitic magnetic forces between the magnetic lay- 
ers can be estimated following Ref.— . The dipole-dipole 
interaction is of the order of attoN at distances of 50 nm 



where i — 1,2 labels each plate, Ui — a is the 
fine structure constant (a — j-), rii = ^ tiiiS) is the 
refractive index of each plate, A = 1 + n| + -I- niX+ 
and: 




It is possible to show that when a = and after a 
little algebra, these coefficients reduce to the ordinary 
Presnel coefficients for a dielectric-dielectric interface^S. 
By its definition, A is always positive and the off-diagonal 
terms have a an overall sign governed by the sign of di. 
Introducing these matrices in ([T]) and assuming that 
both topological insulators are described by an oscilla- 
tor model of the form: 



^(^^)^\^y ^ (9) 

it was showni2, that two different situations could 
arise: 1) when the topological terms had the same signs. 



and the magnetic Casimir forceps is ~ Since we will 

be dealing with distances of the order of \xva. the magnetic 
forces will be even smaller compared to the Casimir force 
that is of the order of pN^^. Thus, we neglect the para- 
sitic magnetic interactions, and focus only on the Casimir 
force experienced by the plates. 



B. Results for isotropic topological insulator plates 

at T = 

In this section we review for completeness the case 
of two isotropic topological insulator plates discussed in 
Refiii^ and extend some of the results. The topological 
part S% in the action does not modify Maxwell's equa- 
tions as long as the constituent relations are changed ac- 
cording to|4]and[5l taken at this stage to be isotropic. 
From these, we can obtain the reflection matrices ^ by 
imposing the continuity of the tangential component of 
H and the normal component of D. For a plate described 
by the optical responses Qi and ei(w) immersed in vacuum 
the reflection matrices are given hySJ^-. 



-1 + n'i + a'i + n^x~ / ' ^ ' 
I 

i.e. sgn(0i) = sgn(02) the Casimir force was attractive 
for all distances, and 2) when the signs were opposite 
sgn(6'i) = — sgn(6'2), it was shown both analytically and 
numerically that a repulsive region at small distances ap- 
peared. At large distances, attraction was recovered and 
the two regimes were separated by a minimum in the 
Casimir energy; an equilibrium point. The dependence 
of the minimum was studied for a one oscillator model, 
showing that the repulsive behavior was favored by re- 
ducing the ratio A detailed version of the proof of 
the existence of a minimum, together with the analytic 
analysis of the large and short distances limits of the 
Casimir energy are given in appendix El 



To end this section we complete the analysis made in 
Ref.^^ by studying the explicit dependence of the posi- 
tion of the equilibrium point {dgq) with 6i and 02- We de- 
fine for convenience the adimensional distance d = —d. 

c 

Whenever the condition sgn(6'i) = — sgn(6'2) is satisfied, 
it is easy to study numerically the dependence of doq, the 
equilibrium position as a function of both 0i 2. The anal- 
ysis is shown in Figl2] In FiglSfa) we plot d^q, against 6*1 
for different values of 62 ranging from tt to Stt. A numer- 
ical fit indicates that for small magnetoelectric couplings 
the dependence is given by: 



(a) 



(b) 



FIG. 2: (color online) (a) d^q as a function of di for different values of 62 ranging from tt to Stt. The solid blue line represents 
deg oc \6id2\-(t>) deq as a function oi 9 = 9i = —62 (black curve with dots). The behavior deviates from 6^ (bottom curve red) 
reasonable for small 6. Beyond 6 = ITvr the equilibrium point is lost (see text). For both figures the remaining parameters of 
the model where fixed to — = 0.45 



4gOc|0i02|. (10) 

For higher magnetoelectric couplings, the behavior de- 
viates from the simple 6'^ \am. This is shovifn in FiglSfb) 
where we plot = 9i = —02- For high values of 9 
the equilibrium position grows faster to larger distances. 
Beyond a critical value, the magnetoelectric part of the 
reflection coefficients overwhelms the ordinary dielectric 
part giving rise to purely repulsive behavior and the equi- 
librium position disappears which is analogous to the 
e(0) — > 1 limit discussed in RefJ^. The inclusion of tem- 
perature will modify further (1101) . as shown in the next 
section. 

As a final comment for this section we mention that the 
obtained dependence of the minimum with 9 suggests 
that an analytic way of deriving this result might be 
possible, especially in the case where 9i = —92- Nev- 
ertheless, the derivation of this result is not transpar- 
ent from this formalism, although other formalisms'^ can 
shed more light to the issue. 

III. CASIMIR-LIFSHITZ FORCE BETWEEN 
ISOTROPIC TOPOLOGICAL INSULATORS AT 

r / 

In this section we include the effect of temperature for 
isotropic plates and we discuss how the repulsive behavior 
can be affected by changes in temperature. 

A. Inclusion of finite temperature effects 

To take into account finite temperature effects within 
the Lifshitz theory one must assume that the system 



is in thermal equilibrium. As a first step, the formal 
replacement^'.^? : 




in ([T|) together with the replacement of ^ with = 
27r-^^Z, the discrete Matsubara frequencies, takes into 
account the effect of temperature. We have defined ks 
to be Boltzmann's constant and the prime denotes that 
the term / = contains a prefactor i compared to the 
other terms in the sum. This set of substitutions imply 
some assumptions which we now briefly discuss. 
First of all we note that formally the dielectric function 
can also have temperature dependence through the pa- 
rameters uJii.,ijJei and 7i defined in Although to 
our knowledge there is no experimental data for topo- 
logical insulators, these parameters are almost temper- 
ature independent for most dielectrics and so we ex- 
clude this effect in our analysis. In addition, it is known 
that all dielectrics have non zero conductivity a{T) at 
T 7^ which modifies the dielectric function through 
e(a;,r) — e{uj) -|-i47r^^. For dielectrics cr(T) depends 
on the band gap A: a{T) ^ exp(— ^^^y ). 

Topological insulators such as Bi2Se3 have A ^ 0.3eV— 
and so we can neglect this contribution as a first approx- 
imation for temperatures below T 3 • 10'^ K which is 
enough for practical purposes. Nevertheless, it is possi- 
ble to show^ that the inclusion of the finite conductivity 
term within the Lifshitz formalism leads to a violation 
of Nernst theorem which states that the entropy should 
tend to zero when T ^ 0. Physically this is related to the 
appearance of a drift current, which leads to Joule heat- 
ing and violates the condition of thermal equilibrium, 
necessary to apply the Lifshitz theory. Hence, the in- 
clusion cr(T) is therefore not justified within the Lifshitz 
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theory and leads to large, unphysical thermal corrections 
which do not account for experimental data^. Either way, 
we are forced to neglect this contribution. Other effects 
of temperature on the coating are discussed in section 

ED 

Having considered these points, we will assume that 
the formal substitution is enough to take into ac- 
count the effect of temperature on the Casimir effect in 
anisotropic topological insulators. 



B. Results at T ^ isotropic plates 

Once established the method to include temperature 
effects we analyze the case of two isotropic topological 
insulating plates described by ([9]) at T 7^ 0. The Lifshitz 
equation ([l} transforms for finite temperatures to: 



1=0 •' ^ > 

(12) 

Where the transverse momentum ^3 is now evaluated 
at discrete frequencies ^; such that = y'fc^ + /c^- It 

is convenient to define the adimensional temperature T: 



T : 



2TrkBT 



(13) 



For frequencies in the optical region, T — 1 correspond 
to temperatures T ~ 2 • 10'^ K, while frequencies in the 
microwave region T — 3.7 correspond to temperatures 
T ~ 300K. Since we will be interested in the latter (the 
topological insulator TlBiSe2 has a resonance frequency^^ 
at ujii = 56cm~^) our plots will range from T = to 
T = 10. We note that the behavior is dominated by low 
frequency resonances as long as higher resonance frequen- 
cies have similar oscillator strengths in © (see auxiliary 
material in Kei^). Because of this, and the absence, 
to our knowledge, of broad frequency experimental data 
regarding the optical response of topological insulators 
(due mainly to their interest as thermoelectrics) we as- 
sume that this is the case and restrict to a one oscilla- 
tor model throughout the paper. Nevertheless, as op- 
tical data becomes available, it is simple to implement 
these techniques and substitute the relevant parameters 
in the model since the overall assumptions will still hold 
(namely, the insulating behavior and the existence of a 
topological magnetoelectric polarizability) . 

We end here the discussion regarding the model and 
the theoretical tools and we now proceed to analyze both 
numerically and analytically the results obtained from 
integrating (H/ 



1. Classical (T ^ 00) and quantum (T 0) limits 

To get some insight into the finite temperature be- 
havior of the system one can ask what happens to the 
Casimir force in the (formal) limits when T — 00 and 
T 0. Trivially, if one first sets T — the expres- 
sion ([T|) is recovered and the analysis of Ref.^^ still holds. 
Hence a region of repulsion is expected to exist if the 
topological magnetoelectric polarizabilities of the plates 
have different signs. For high temperature, T — > 00 also 
referred to as the classical limit, all the terms in the / sum 
above are exponentially suppressed except for the I = 
terni^. Thus, for isotropic plates, the dielectric function 
can be replaced by the zero frequency dielectric function. 
Let us derive the reflection coefficients in this case. 
Directly from ([7]), and taking into account that nf = 
the reflection coefficients at zero frequency take the form: 
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1 

2a 



(l-e(0)-a2 + y^X_(0,fc||)) (14) 
(-1 + e(0) + 52 + V<0)x- (0, fc||)) (15) 

(16) 



where A;=o = 1 + e(0) + + ye(0)x+(0, fc|| ). From 
their definition ([5]), the functions x± zero frequency 
can be expressed in terms of e(0) only (fc|| dependence 
cancels in this limit): 



1 



e(0)' 



(17) 



Bringing together all simplifications the reflection co- 
efficients ([7]) take the form: 



R, 



2a 



2(e(0) + 1) + a2 V 2a 2(e(0) - 1) + a^ 



(18) 



We emphasize that although similar in structure, these 
reflection coefficients are different to the ones obtained at 
zero temperature and d — >■ c» in appendix VK\ bv rescaling 
the Casimir energy density with d. 

One can introduce this expression in (jl2p and study the 

behavior of the Casimir energy as a function of the two 

2 

parameters of the model: e(0) = 1 + ^ and Oi_2- In this 
case one can rescale the integration variables with d and 
focus on the force defined as F = —ddEc{d) which takes 
the form: 



Fciid) 



d^-K 



f{e{0), 61,02), 



(19) 



where / is a complicated integral expression of the pa- 
rameters (see appendix |B] for details). It is evident from 
the functional form of (IT51) with d that in the classical 
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FIG. 3: (color online) Attraction versus repulsion in the 
classical limit (T — > oo) as a function of — and - where 
9i = —02 = d. The white upper region shows attraction while 
the lower green region shows repulsion. The black region is a 
forbidden region for the parameters due to the positive energy 
condition < \J e(0) (see text). 



limit the Casimir force is either attractive or repulsive 
for all distances, depending only on the sign of /. Thus, 
in contrast to the zero temperature case, no equilibrium 
point is expected. 

When sgn(6'i) = sgn(6'2) no repulsion is obtained for any 
values of the parameters, as expected from the results 
in Reffi^. Whenever sgn(0i) ^ sgn(6'2) it is possible to 
draw the diagram of Fig. |3] where the attractive and 
repulsive behaviors are shown as a function of the two 
parameters for the particular case where 9i — —92 = 0. 
To understand this diagram one first has to note that 
the topological magnetoelectric polarizability has an up- 
per bound to satisfy the positive energy conditiort^ for 
magnetoelectric materials. In the present context, this 
condition implies that < -\/e(0) which excludes the 
black region on the right of Fig. [3l The other two regions 
indicate attractive (white top region) or repulsive (green 
low region) behavior. From this diagram it is transparent 
that even in the classical limit of very high temperature 
there exists a region of parameters where complete re- 
pulsion is obtained. It is relevant to point out that a low 
dielectric response is needed for this to happen. High 
topological magnetoelectric polarizability also enhances 
the repulsive behavior. 
Can one understand the origin of repulsion in the classi- 
cal limit analytically? The answer is in fact affirmative. 
For small 9i,2 and 1 — e(0), it is possible to write a closed 
analytical form for the Casimir force of the parameters 
although for the sake of clarity we leave this to appendix 
IB] In addition to this, we have included a discussion of 
the classical limit in terms of the Casimir energy den- 
sity (see appendix [B]) in the spirit of reference^^ study- 
ing the relative strength of the diagonal and off-diagonal 



reflection coefficients. Both approaches lead to the same 
result: repulsion is possible in the classical limit for a cer- 
tain range of parameters which are summarized in Fig|3l 
Having studied these limiting cases we know proceed to 
the general case where T ^ 0. 

2. T 7^ 0; the general case 

The effect of increasing temperature for a situation 
where sgn(0i) = — sgn(6'2) is presented in FigH) In gen- 
eral, temperature acts against repulsive behavior, driv- 
ing it towards shorter distances as temperature increases. 
The evolution of the equilibrium position with respect to 
the adimensional temperature T is shown in Fig. lUJb). 
For high values of T (T > 5) it is possible to complete 
equation (fTO|) (vahd for small values of ^1,2) to write: 



This simple relation, extracted from numerical data, 
shows that for small topological magnetoelectric charge 
and high temperature both effects compete as regards to 
the equilibrium position. 

One can summarize the results for this section as follows. 
In the classical limit, one can have either attraction or 
repulsion, regardless of the distance. Starting with at- 
traction at the classical limit, i.e T — > 00, decreasing 
temperature will develop a deeper minimum which will 
shift to longer distances as shown in FigHJa) and (b). 
Consistent with the discussion presented the end of sec- 
tioniniB, in some special cases where the magnetoelectric 
is bigger than a critical value one can end with repulsion 
at T = without a minimum. 

Starting from repulsion in the classical limit will not de- 
velop a minimum when decreasing T but it will shift the 
curve of repulsive behavior towards larger distances to 
end with repulsion in the quantum limit. This is shown 
in FigUfc). It must be noted that for this to happen, a 
very low e(0) is needed when ~ tt and as 9 is increased, 
higher e(0) are capable of giving rise to the complete re- 
pulsive behavior. Therefore we expect this not to be the 
relevant situation for real topological insulators, which 
have low 9, but rather the one presented in Fig. |H[a). 
An extension of this analysis could be relevant for other 
magnetoelectric materials such as Cr203 that can present 
higher axion coupling as well as more general magneto- 
electric couplings'^^ . 

We note for completeness that when both topological 
magnetoelectric couplings have equal signs, temperature 
acts to increase attraction as expected from ordinary di- 
electrics. 

We finish this section with a technical comment regard- 
ing the commutativity of the limits. We have discussed 
the high and low temperature limits, and studied the dis- 
tance dependence numerically. At this point, one might 
ask whether the limits commute, i.e. if one evaluates first 



'''' (b) (c) 

FIG. 4: (color online) (a) Casimir energy density (in units of Eo — Ahc/{2Tx)^{u)B,/cf') as a function of the dimensionless 
distance d for ^ = 0.45 for temperatures ranging from T = — 10 and 9i = —62 = Stt. Temperature pushes the repulsive 
behavior to shorter distances, (b) Equilibrium distance as a function of temperature, (c) Casimir energy density (in units of 
as a function of d for ^ = 0.16 for temperatures ranging from T = — 10. For values where there 
is repulsion, decreasing the temperature shifts the curve towards smaller distances as shown inside the inset figure. 



the (i -> or d — !> 00 limit first and then consider T ^> 
and T — i> 00, would the result remain as stated above? 
The fact is that in some cases, the limits commute while 
in others they do not. In particular, is clear from ex- 
pression (IT^ that the limits d — J> 00 and T 00 will 
commute. It is easy to show that, in either order, both 
result into the same reflection coefBcients given by (IT51) . 
For the case when d — and T — > a similar situation 
occurs and the limits also commute. The interesting case 
is however when d — > and T ^ 00 (or viceversa). If 
one expands the exponential function, then taking the 
temperature to infinity breaks down the expansion and 
this limit is not reconciled with the T — 00, d — > which 
typically gives attraction. Although it might seem sur- 
prising at first, one should recall other physical situations 
where a similar situation happens, one particularly im- 
pressive example taking place when calculating the DC 
or minimal conductivity of graphene where the order of 
evaluating the limits is crucial (for a recent review of the 
problem see Ref.'^^ section IV. A). Nevertheless, it should 
be noted that, with the numerical results in mind, there 
is a way around the puzzle by imposing that we take the 
limits in such a way that the product Td is a constant. 
In this way, the / = is the leading term in the sum (fT2|) 
and we recover the reflection coefficients (fT8|) relevant for 
high temperature, making the limits commute. 



IV. CASIMIR FORCE BETWEEN 
TOPOLOGICAL INSULATORS WITH UNIAXIAL 
ANISOTROPY AT T = 

In this section we again set the temperature to zero 
and focus on the effect that uniaxial anisotropy has on 
the repulsive behavior. 



A. Anisotropy in topological insulators 

For a given material its linear response to an electro- 
magnetic wave, or in turn, its reflection coefficients de- 
fined in ([1]) depend strongly on the dielectric tensor e,;j, 
the magnetic susceptibility tensor /i^ and the magneto- 
electric tensor 9ij. Since we are dealing with non mag- 
netic topological insulators we will assume that fiij — Sij , 
although we shall derive the reflection coefficients for the 
general case in which the susceptibility is considered. 

Crystal symmetry defines the other two tensors. Due 
to time reversal symmetry, for topological insulators the 
magnetoelectric tensor can only have the form O^j = 9Sij 
with 9 — TT (mod(27r)) for topological insulators and 
9 = for ordinary insulators^. It is interesting to note 
that for general magnetoelectric crystals, other magne- 
toelectric couplings can be present in the action ([3|) , and 
9ij could have a more elaborate form'^^'^^ which could in 
principle affect the sign of the Casimir force. We will re- 
strict ourselves to the case of a constant isotropic magne- 
toelectric tensor 9ij = 9Sij, relevant for crystalline topo- 
logical insulators^. 

The form of Cij is also determined by crystal symmetry 
and one can distinguish the cubic, uniaxial and biaxial 
crystals^. The simplest case, the cubic structure, has 
eij = e(u!)Sij, which is just the isotropic case consid- 
ered in previous sections. Although interesting as a first 
approximatioE^i^, this model is incomplete for topologi- 
cal insulators. The prototypical examples of topological 
insulating crystals BiaSes, BizTea , SbaTea^iiSe^ 
novel Thallium based III-V-VI2 compound o^^'^^" — have 
rhombohedral symmetry, which fall into the category of 
uniaxial crystals and hence we will focus on this case, 
leaving biaxial crystals aside. Uniaxial crystals have a 
definite optical axis which coincides with the principal 
crystal axis. The dielectric tensor can thus be written in 
the form: e^j = diag(e^, ex, £2), where the optical axis 
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is aligned with the z axis. In ah our calculations, the 
optical axis of the two Casimir plates are supposed to be 
aligned so that no torque is induced in the system^. 



B. Fresnel coefficients for topological insulators 
with uniaxial anisotropy at T = 

Fresnel coefhcients are derived from Maxwell's equa- 
tions by imposing continuity relations for the normal 



J 



components of D and tangential components of H as 
discussed above for the isotropic case (see appendix [C|. 
After some tedious but straightforward work one ob- 
tains the reflection matrices for an anisotropic topological 
insulator-vacuum interface: 



R, 



2aq"fij_k^ ie±k^ 



2aq" j.i±kz 
q"){^i±k, + q')+q"k,^i^a^ 



(21) 



where Aa„ = {fi±kz + q') {e±kz + q") + q"kzfi±a'^ , 

, ,2 , ,2 ,, , ,2 



^1 - ~^ ^jj, 9 



/2 _ ^ , ,, l2A^-L 1,2 _ ~ ^ 



1? 



fcy — and a — aO/i:. To compute the Casimir energy, 

one should define these in the imaginary frequency axis 
and perform the substitution lo = i^. We stress that 
these reduce to the ordinary anisotropic coefficients pre- 
sented in appendix [Cl when a = 0. They also reduce to 
the isotropic reflection coefficients (O when /i and e are 
isotropic. 

To investigate the effect of anisotropy it is necessary to 
model the dielectric response of the material. In general, 
insulators can be described by its resonance frequencies, 
hence we model e± and with a sum of oscillators: 



i 



(22) 



(23) 



We have evaluated the dielectric function at imaginary 
frequencies, as demanded by equation ([T]). The parame- 
ters ujej^ i and ; indicate the strength of the oscillator 
while UJR^ . and ujn^ . account for the resonance frequen- 
cies of the oscillators. Finally, we include the damping 
parameters 7^,^ and ^jz^i for completeness although they 
do not play a mayor role on Casimir physics. 

A multiple oscillator model does not significantly alter 
the general behavior of the Casimir force^^ and so, in 
order to clarify the effects that different parameters have 
on the Casimir energy density we will restrict the models 
([2^ and to one oscillator models. In this scheme we 
will investigate two aspects: (i) the effect of modifying 
the relative strength of the parallel component against 
the perpendicular i.e. cje, over ojg^ and (ii) the effect of 
the relative movement of the position of the resonance 
frequencies ujjf^ and ujji^ . 



C. Effect of the relative strength of the parallel 
component against the perpendicular component 

To study this effect we simplify equations ([^ and (^5)) 
to unmask the effect under study (i.e. point (i) above). 
Redefining the parameters in units of lor^ : 



ezm = 1 + 



2 



(24) 
(25) 



In this case we have chosen ujn^ = ujr, and 7 = — 
7_R^ = 0.01 in units of ujr^ . To study the effect of chang- 
ing the relative strength of the dielectric components we 
may fix = 0.45 (in units of wjjj^) and change Wgj^ in 
the interval — 1. The results are shown in FiglS] We 
have fixed 9i — —62 = (i.e. the topologically trivial 
anisotropic case) in Fig. [Sfa), 9i — 62 = tt in Figure [5jb) 
and 9i = —62 = tt in Fig. [Sfc). 

Fig. EJa) is a reference figure where only anisotropy 
is studied with = 0. In this case, attraction occurs 
since 6i = Q and there is no mixing of polarizations since 
the diagonal component in (PT|) is zero (see appendix 
As cje^ grows the tendency is to favor attraction and to 
increase the absolute value of the Casimir energy at a 
given distance. In Fig. [S^b) the magnetoelectric terms 
have equal sign and so no repulsion is obtained just as 
it is expected from the isotropic casein. When 0i and 02 
have opposite sing, Oi = —02 = tt, the Casimir energy 
develops a minimum, shown in Fig. [5jc) . The actual dis- 
tance at which the minimum appears is shifted to shorter 
distances as oJe^ is increased. This enables to draw the 
conclusion that for the minimum to shift to larger dis- 
tances, or in other words, for the repulsive behavior to 
be enhanced, the material in question has to favor the 
strength of oscillators in the direction parallel to the op- 
tical axis. 
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(a) (b) (c) 

FIG. 5: (color online)Efrect of changing the relative oscillator strengths between parallel and perpendicular components of the 
dielectric function on the Casimir energy density (in units of Eq = Ahc/ {2tv)'^ {ujr^ /c)^) as a function of the dimensionless 
distance d with (a) 61,2 — and (b) 61 = 62 = tt (c) 61 = —62 — tt. (a) Attraction is obtained as no magnetoelectric term 
is included, (b) Both magnetoelectric couplings have the same sign and so no repulsion occurs, (c) Repulsion appears at 
short distances. Increasing oje^ makes the minimum shift to lower distances while increasing its depth. Therefore, repulsion is 
favored by the condition lj^^ > uj^^ 



Physically, repulsion in this system is due to the mixing 
of polarizations^^. Thus, suppression of repulsion can be 
traced back to the relative suppression of the off-diagonal 
terms in (PT|) versus the diagonal terms just as in the 
anisotropic case studied in^^. When increasing e± fixing 
ez to a constant, the off-diagonal terms vanish, confirm- 
ing the numerical analysis. 



D. Effect of the relative position of the oscillator 
frequencies 

To study this case, the dielectric functions of both com- 
ponents can be written as: 



1 



(26) 
(27) 



In this case we have chosen — t^e^ = — 0.45 
and as before 7 = 0.01 (in units of Wfl^). The param- 

eter fi = -^^^ determines the relative position between 

both oscillators: for /3 > 1, has a higher resonance 
frequency than and for /3 < 1 the opposite situation 
occurs. Again we have fixed 9i — —62 = in Fig. El^a), 
6*1 = 6*2 = = TT in Fig. mb) and 9i = -62 = 9^ = tt in 
Fig. He). 

As in the previous case, Fig. [D^a) a reference figure 
where only parameter /? is modified. We have chosen 
the values of /? ranging from /? — 10^^ to /3 = 10'^, so 
a broad range of situations are represented. In this sit- 
uation, attraction occurs since 9i = and there is no 
mixing of polarizations. As /3 grows the tendency is to 
favor attraction and to increase the absolute value of the 



Casimir energy at a given distance. In Fig. the mag- 
netoelectric terms have equal sign and so no repulsion is 
obtained. Fig. [6jc) shows how the minimum is shifted to 
lower distances as lurj^ is increased, i.e. /? is increased. 
Hence, for the repulsive behavior to be favored, the ma- 
terial in question has to favor low resonance frequencies 
in the parallel direction, while high frequencies in the 
perpendicular direction. This is consistent with the be- 
havior discussed above where suppression of ej_ enhances 
repulsion. In this case, increasing /3 suppresses e± and 
favors repulsion^. 

We summarize the results of this section in what follows. 
From the numerical results one can readily infer what 
conditions are necessary for repulsion to be observed, 
given that the topological magnetoelectric terms of the 
two materials have opposing signs. For the repulsive be- 
havior to appear at the largest distances possible, one 
should search for a material where the direction parallel 
to the optical axis has a bigger oscillator strength than 
the corresponding strength of the oscillator in the par- 
allel direction cje, > Wgj^ . Further enhancement of the 
repulsive behavior can be achieved if the oscillator reso- 
nances in the direction perpendicular to the optical axis 
are at higher frequency than resonances for the parallel 
direction, i.e ujn^ > ujr^ . In both situations, ej^ is sup- 
pressed relative to and the mixing of polarizations is 
maximized. 

With these numerical calculations we have established 
the main directions to enhance repulsive behavior with 
topological insulators at low temperature. We now dis- 
cuss the effect of temperature on these findings. 
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(a) 



(b) 



(c) 



FIG. 6: (color online) Effect of changing the relative positions of the oscillator resonances between parallel and perpendicular 
components of the dielectric function on the Casimir energy density (in units of -Bo = AHc/ {27r)'^ {bJn^ /c)^) as a function of 
the dimensionless distance d with (a) Si, 2 = and (b) Oi = 62 = n (c) Oi = — 6*2 = vr. (a) Attraction is obtained as no 
magnetoelectric term is included, (b) Both magnetoelectric couplings have the same sign and so no repulsion occurs, (c) 
Increasing lor^ makes the minimum shift to higher distances while decreasing its depth. Therefore, repulsion is favored by 
increasing /3.(see text) 



V. CASIMIR FORCE BETWEEN 
TOPOLOGICAL INSULATORS WITH UNIAXIAL 
ANISOTROPY AT T 7^ 

The result of including the two effects subject of this 
work, temperature and anisotropy together is shown in 
Fig [71 To illustrate that the repulsive behavior is still 
present we study two different cases. In the first of them, 
corresponding to Fig. [TJa) we vary the adimensional tem- 
perature from T = to T = 7 while fixing the uniaxial 
anisotropy parameters and 61 = —62 = tt. In particular 
we fix the oscillators strengths to lu^^ — 0.45, ojej^ — 0.3 
and the relative position of the resonances to /3 = 5, 
i.e. the perpendicular oscillator is shifted by a factor of 
five with respect to ajj?^ . As expected by previous sec- 
tions, the figure proves that temperature works against 
the repulsive behavior, shifting repulsion to shorter dis- 
tances as temperature is increased. In a second situation 
we illustrate the effect of varying anisotropy at a given 
temperature. In this case, we fix the temperature to 
T = 5, the oscillator strength in the parallel direction to 
= 0.45 and /3 = 5 while varying lj^^ from uj^^^ = to 
= 1- The results show that enhancing perpendicular 
response, even at finite temperature, favors repulsive be- 
havior. Consistently with previous sections, equal signs 
of 9i^2 the system returns to attraction at all distances. 



VI. DISCUSSION AND CONCLUSIONS 

In this work we have studied the effect of anisotropy 
and temperature on the repulsive Casimir effect with 
three dimensional topological insulators. We have found 
that the repulsive behavior is still present even when in- 
cluding these effects. In particular, increasing temper- 
ature tends to reduce repulsion, similar to what is ex- 



pected from different situations such as metamaterials^. 
On the other hand, enhancing the optical response in the 
direction parallel to the optical axis of the topological in- 
sulator works in favor of repulsion. This can be achieved 
either by increasing the relative strength of the oscillator 
in the parallel direction or by searching for resonance fre- 
quencies which suppress the perpendicular response. We 
have also determined that both effects together still allow 
for repulsion. All these effects can be understood as an 
enhancement of the off-diagonal terms in the reflection 
coefficients, which favors polarization mixing and thus 
repulsion, whenever the signs of the topological magne- 
toelectric terms are opposite in both plates. 
Although these specific conclusions are reached, it is dif- 
ficult to draw conclusions for real materials since the op- 
tical response of topological insulators still needs further 
experimental characterization. In particular, values for 
the optical parameters are still lacking, in contrast with 
other materials which are traditionally useful for measur- 
ing the Casimir force^i^. Nevertheless, in this work, we 
have aimed also to establish the theoretical framework for 
handling anisotropy in topological insulators in the con- 
text of their optical properties, which to our knowledge 
was absent in the literature. The reflection coefficients 
for anisotropic materials can also be helpful in achieving 
the optical characterization of these materials and could 
be relevant for metrology purposes, as suggested in pre- 
vious works'*^'—. This theoretical work on anisotropy 
might also be used as a starting point for the study of 
the Casimir force in situations involving other magneto- 
electric materials. 

A flnal issue left out in this work is the effect of temper- 
ature on the magnetic coating. Even though magnetic 
interactions where shown to be irrelevant temperature 
will have an effect on the magnetization. In particular, 
the magnetic moments will fluctuate with increasing tem- 
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(a) (b) 

FIG. 7: (color online) Casimir energy density (in units of Eo = Afic/ {2-11)^ [ojr^z / cf') and 9i — —O2 — vr including anisotropy 
and temperature effects: (a) Effect of changing the temperature from T — — 7 with ijje,z ~ 0.45, cje,± ~ 0.3 and /3 = 5. (b) 
Effect of varying uje,i_ =0 — 1 with T — 5, uje.,z ~ 0.45 and /? = 5. 



perature and destroy ferromagnetism when the system is 
over the Curie temperature (Tc) of the coating. This tem- 
perature depends on the specific coating material which 
should be an insulating or semiconducting ferromagnet 
with as high as possible Curie temperature. Although 
ferromagnets are usually metallic, examples of insulating 
ferromagnets exist, for example oxides such as the rare 
earth oxide EuO with Curie temperature close to 70K— , 
and have been synthesized and used in the past, such 
as the ferromagnetic insulator GdN. The latter has been 
shown to have a Curie temperature close to 60K"'^, it can 
be grown as a thin film of thickness ^ Gnm"*^ and could 
be a suitable candidate for the ferromagnetic covering. 
In conclusion we expect that the tools and behaviors 
studied in this work might lead towards a complete opti- 
cal characterization of topological insulators, a first step 
towards Casimir like experiments in which the descrip- 
tion of factors such as temperature and anisotropy should 
be relevant. 



Appendix A: Proof of the existence of the minimum 
at T = 0: Large and short distance limits of the 
Casimir Energy density: 

In this appendix we give detailed information and ex- 
pressions of the some results presented in^'^ so this work 
is self-contained. We prove the existence of the mini- 
mum by discussing the high and low frequency limits of 
the Casimir energy density. The Casimir energy stored 
between the plates is given by: 



Ecjd) 
Ah 



1-K 



d^kll 



logdet [1 - Ri • Rae^^fcsdj 



(Al) 

where Ri is the 2x2 matrix defined in ([7|). To prove the 
existence of the minimum we assume that the refiection 
coefficients are given by the expression ([7]) and that the 
dielectric function satisfies the following two analytical 
properties: 1) finite dielectric permittivity at zero fre- 
quency (e(0) < oo) and 2) high frequency transparency: 
e(a;) — 1 when a; — > cxd. For an insulator one can assume 
the dielectric function is given in (jO]). For simplicity let's 
assume that: 
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(A2) 



although the derivation does not depend on the ana- 
lytical form of the dielectric function as long as it fulfills 
the mentioned conditions (the proof holds for more than 
one oscillator). For analytical traceability we assume the 
particular situation where 9\ = — ^2 (where labels 1 and 
2 identify the Casimir plates). 

The first step is to rescale Ec{d) with d. If in the expres- 
sion for Ec{d), ^ and kii are rescaled to contain d, tan 
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overall factor 1/d? appears in front of the integral. The 
reflection matrices are to be evaluated now at the rescaled 
frequency and momenta £,/d and kii/d. Thus, ah of the d 



J 



dependence can be transferred to the reflection matrices 
and the overall prefactor. The expression for the 
energy is: 



Ea{d) 
Eo 



I d^ I d^k||logdet 



1 -Ri 



r 



d' d 



R2 



d 



-2k3 



(A3) 



where fcg = ^ + fcjj is defined through the rescaled 
variables (which include c the speed of light in vacuum 
and lob) and Eq — Ahc/ {2tt)'^ {w r/ c)^ where A is the area 
of the plates. The key point is that reflection matrices 
are evaluated at these rescaled variables, with a rescaled 
dielectric function: 



,{t^/d) = 1 + 



1 



ijjji d 



(A4) 



The two conditions imposed to e(i^) ensure that the 
reflection matrices are not singular in the limits ^/d ^ 
and ^/d 00. Hence, Ec{d ^ 0) ^ ±00 and Ec{d — >■ 
00) ^ 0. 



J 



The way the integrand approaches these limits deter- 
mines the sign of Ecid). For instance, if the integrand 
is positive at small distances and negative at large dis- 
tances, necessarily a minimum exists at an intermediate 
distance. When 61 — —62, this is exactly what happens 
(unless e(0) = 1, where both limits are positive and hence 
long range repulsion is obtained), as it will be shown in 
what follows. 

We now evaluate the integrand in (|A1[) . When 9i = 
—92 = 9 the reflection matrices describing both topolog- 
ical insulators can be written as: 



R 



± — 



±rsp(j^,k||) rp(i^,k||) 



(A5) 



Introducing this inside (|A1[) the integrand follows: 



/ — log det 



1 R I R 



-2k\ 



log 



-2k\ 



r 



+ e 



-4k\ 



sp 



Notice that the last term although always positive, will 
play no role in what follows since it is always suppressed 



over the first term (note that 



sj ' pi ' sp 



<1). 



cies) we have transparency, n^i^/d) = e(^/(i) ^ 1 then 
we see that x_ and: 



The limit of short distances: 

For d ^ and using the high frequency transparency 
of the dielectric function we now show that \ri\ << \rsp\ 
{i = s,p): Notice first that the denominator A defined in 
([7]) is common to all terms so it cannot play a role on the 
relative magnitude of the coefficients. We hence study 
the behavior of x_ at small distances, given by: 



2 + +X+ 



and 



2a 



2 + a^ +X+' 



d' d 



e 



Remembering that at small distances (large frequen- 



(A6) 



since the first are of order O(a^) and the second are 
of order 0{a) (remember that a is proportional to the 
fine structure constant a) we have that \ri\ << \rsp\ (i — 
s,p). Hence the integrand is positive (since the integrand 
has the form / = ln(l + A) where A > 0) and so Ec{d 
0) +00. 
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The limit of large distances: 

When d — > oo the reflection coefficients take the form 

1 - e(0) - a2 + v/^X- 



for the diagonal part (with a similar expression for rp) 
and 



2a 



for the off diagonal. We have defined the quantity 
e(0) = 1 + (^-^^ ■ In this case, depending on the values 

of e(0) different behaviors emerge. Since e(0) > 1 we can 
distinguish to extreme limits, one where e(0) = 1 and the 
other with e(0) >> 1. 

When e(0) = 1 one can see that we return to the pre- 
vious case since the quantity X- in this limit also goes 
to zero. Hence \rsp\ » \ri\ [i — s,p) is satisfied for all 
distances and Ec{d) is always positive .Therefore, using 
the fact that that Ec{d) when d ^ oo and that 

I 



Ec{d) — J> +00 when d — ?• we deduce that there is no 
minimum and that the force is always repulsive. This 
was confirmed by the numerical calculations presented 



When e(0) >> 1, we see that the opposite condition, 
\fs.,p\ >> {fspl, is satisfied. Even in the worst case when 
X- is smallest, e(^) in rs is always larger than 2a in rsp. 
The integrand at large distances is a negative quantity 
(since the integrand now has the form / = ln(l — B) with 
B > 0) and so Ec{d) approaches to zero from negative 
values {Ec{d) —oo for d -> cxd). Since at small dis- 
tances Ec{d) — > +00 we conclude that there must be a 
minimum at an intermediate distance < dm < oo, since 
the function must cross the d axis. 

Notice that when e(0) is strictly infinity we recover the 
case of an ideal metal where r^.p = Tl a-nd rsp = 0, with 
attraction at all distances. 

To sum up, as we increase e(0) from one, a minimum 
develops at an intermediate distance dm- This distance 
shifts to lower values as we increase e(0) until, at e = oo 
we recover the ideal metal case where complete attrac- 
tion occurs. 

In the case where 9i — 62 the signs inside / change to 
give: 



I = log 



{2r 



The predominant term inside the logarithm is always neg- 
ative and hence the integrand is always negative. There- 
fore at small distances Ec{d) — J> —00 and at large dis- 
tances Ec{d) — >■ approaching this limit from nega- 
tive values, recovering attraction at all intermediate dis- 
tances. 

Finally, it can be checked by analogous methods, that 
the case where one 9 is zero and the other one is not 
(dielectric - topological insulator case) results in Casimir 



J 



attraction for all distances. 



Appendix B: Classical limit: analytical expressions 

The Casimir force at non zero temperature is defined 
as F = -ddEaid) with Ec{d) given by ([72]): 



F(T,d) 



dk\ 



fc||/c3Tr 



1 Ri R2e 



Ri Rae 



(Bl) 



r 



where ^3 = ^/kf. +T'^P and Ri are the reflection coef- 



.•3-y..|| 

ficients of each plate. In this expression the momentum 
variable is rescaled in units of and the integral de- 
pends only on d and T. As discussed in the main text, at 
high temperature, only the I — term contributes since 



the rest are exponentially suppressed. Keeping only this 
term in the integral one can recast this expression as: 



Fciid) 



d'^TT 



(B2) 
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Integrating this expression for different values of 9i^2 
and e(0) one obtains the diagram in FiglJ] (whenever 0i 2 
have opposite signs but equal magnitude 9). 
In order to derive analytical results further approxima- 
tions must be considered. In the rest of this appendix 
we will assume that we are always in the region of small 



magnetoelectric coupling a << 1, which is justified as 
long this coefficient is of the order of the fine structure 
constant, assumption to be expected in real topological 
insulators. Hence, for small topological magnetoelectric 
couplings we can evaluate the integral to give: 



Fciid) = - 



Li?. 



m + ly 



a2 + 2aia2 



6(0) + 1 



6(0) 
I 



1 



Ij (6(0)-l)(6(0) + l) 



i«2 



fm- 









(B3) 



where a = ^ and Li„ is the polylogarithm of order n. 
When e(0) — ^ 1, the force takes the simple form: 



Fci{d) 



Wd^TT 



aia2- 



(B4) 



From this last expression it is transparent that one can 
still get repulsion in the classical limit when the signs of 
01 and 02 are opposite. 

As anticipated in the main text one can analyze the 
Casimir energy density in the classical limit in the spirit 
of Refii^, that is, exploring the relative magnitudes of 
the diagonal and off-diagonal parts of the reflection co- 
efficients. In this case, we assume for simplicity that 
both topological magnetoelectric polarizabilities have the 
same magnitude but can have equal or opposite signs, i.e. 
sgn(0i) = ±sgn(02). If we introduce these reflexion ma- 
trices in expression ()12|) and take the high temperature 
limit (l — 0), after some algebra one arrives to: 



sgn(6'i) = - 
makes E cc 



sgn(^2) and 6(0) ^l,a;<0so/>0 which 
and the force is repulsive since a nega- 
tive (positive) force, or equivalently a positive (negative) 
slope of Ec{d), corresponds to attraction (repulsion) of 
the plates. 

Indeed, when one neglects the second term inside the log- 
arithm in (jBSP the integral can be expressed in terms of 
a Polylogratihm: 



ksT 
167rd2 



Li3(A±) 



(B6) 



The function 'L\?,{A±) has the sign of A±. Therefore, 
when A± > 0(< 0) we have attraction (repulsion). Using 
the fact that Lis (a;) x + 0{x^) when x ^ and that 

^ when e(0) — 1 and to leading order in 
we can readily obtain the Casimir force from (|B6p which 
gives: 



ksT 
And^ 



■ / log [1 - A±e-2'=ii + Be-^^ii] , 

Jq 

(B5) 

where the ± corresponds to each of the two cases 
sgn((?i) = ±sgn(6'2) and we have rescaled the integral 
variables with the distance. The functions A± and B are 
defined as: 



B 



D 
D 



(4(6(0) - 1)2 + (4(6(0) - 1) ± 8)a^ + 2a^), 



where D = (2(e(0) -I- 1) -I- d^) . Note that these two 
coefficients depend on the zero frequency dielectric func- 
tion 6(0) and the absolute value of the topological mag- 
netoelectric polarizability \9\ as well as its sign in the 
case of A± . These functions govern the sign of the inte- 
gral. Function B plays a secondary role since it is expo- 
nentially suppressed. It is clear, that when 6(0) — J> 1 
the integrand takes the form / = ln(l — x) where x 
is positive or negative depending on A±. For A_, i.e. 



Fci{d) 



keT . 
led^Tr' 



(B7) 



This equation is exactly equation (jB4[) when di = 
—d2 = OL. Consistently, both approaches are equivalent 
and confirm our numerical results, where we found that 
even at high temperature repulsion can occur. Neverthe- 
less, for more realistic situations where higher values of 
6(0) are expected, there is a competition between param- 
eters, which is represented in FiglS] Only high values of 
can compete with large enough 6(0) and thus we expect 
attraction to occur at high temperatures. 



Appendix C: Fresnel coefficients for topological 
insulators with uniaxial anisotropy 

In this appendix, the Fresnel coefficients for topologi- 
cal insulators with uniaxial anisotropy will be discussed. 
First, we will briefly review the derivation of the reflec- 
tion coefficients for an ordinary dielectric - vacuum in- 
terface so that the extension to the topological insulator 
- vacuum interface is more transparent. 
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Presnel coefficients for uniaxial material plates 

As an illustrative case, before adding the axionic term, 
we proceed to solve Maxwell's equations for an uniax- 
ial material which will demonstrate the procedure to 
follow in the next section. We start with an uniax- 
ial material with its optical axis parallel to the plate's 
normal. The corresponding dielectric permittivity ten- 
sor is given by eij — diag(e_L, e_L, e^) and, for complete- 
ness, the magnetic susceptibility tensor will be included 
as = diag(/i±,/ix,^z). The subindex z indicates 
that the optical axis is aligned with the z axis, chosen 
also to be the surface normal. The procedure is to solve 
Maxwell's equation in vacuum, then inside the uniaxial 
media, and finally impose boundary conditions to deter- 
mine the reflection amplitudes. The first step is therefore 
to solve Maxwell's equations in vacuum by proposing a 
plane wave solution of the form (sec for instance Refi^) 



c 



For the incoming wave, and 



E7. 



R\\y + — (fcxX + k^z) ) e 

id 



i(k^x — kzZ—ujt) 



for the reflected wave, where we have used that fc™ = 
—k^^-^ with obvious notation. The cartesian unit vectors 
are represented by x, y and z and we define k^ — ^sin^j 
and kz = ^cosOi, where 6i is the angle of incidence. 
The problem consists in finding the relative amplitudes 
and R±,R^^. Their quotients will define the en- 
tries of the reflection matrices ([2]) . 

The second step is to solve Maxwell's equations inside 
the uniaxial material. Since there is translational invari- 
ance along the interface, kx must be conserved so the 
transmitted wave can have the form: 



= Q^i{qz+k^x~ujt) ^ 
JJ^ _ Y^^i{qz+k^x-u>t) 



(CI) 
(C2) 



where q is the transverse transmitted momentum to be 
determined by finding the dispersion relation. From the 

1 aB 

Maxwell equation V x E = with Bi — ^J-ijHj 

one can obtain the conditions for vectors e(z) and h(z) 
which are: 



ck 

hz = e 

cq 



UJfJ.± 

{kxCz - qex) ■ 



(C3) 



From the Maxwell equation V x H = — — with Di 

c at 

e-ijEj one obtains: 



ckx , 

e -^h 

UJ€± 



(C4) 



From the sets (|C4p and (|C3[) , after a few steps one can 
see that if Cy ^ 0, ^ and hy — Cx ^ 0, q must 
satisfy: 



V C / flz 



(C5) 



This is just the transverse electric mode inside the ma- 
terial. Then, in this case: 



ey - tlx 



(C6) 



In a similar fashion if Cy — 0, hx — and hy =i 
ex ^0, q must satisfy: 



q^^q"^=(^)\^e^-kf-^, (C7) 



and: 



^ - cq" 



(C8) 



which corresponds to the transverse magnetic mode in- 
side the uniaxial crystal. 

It is now time to impose the boundary conditions. As 
mentioned above, tangential components of E and H 
must be continuous along the interface, thus: 



(^11 -R\\)- = ex, 
(jj 

{A± + R±) = ey, 

(A|| - hy, 



(i?± - 



(C9) 

(CIO) 
(Cll) 
(C12) 



The last two equations can be further simplified by 
using dCHI and jCHI) to give: 



(A||+i?|| 
(R^-A^)- 



cq" cq" 



^11)-' 



cq' 
( 



cq' 



(A^-R^). 



From these it is clear that the transverse electric and 
the transverse magnetic mode are decoupled, a well 
known result from ordinary electromagnetic theory^S. 
One can now solve for the quotients R±/A± and ^||/^|| 
and find the ordinary reflection coefficients for a uniaxial 
crystal, which in matrix form define the refiection matrix: 



16 



(r) 

s 

(r) 




^e_L/i_L - fcy 



(i) 

s 

(i) 



(C13) 




where Es.p and i?s_p are reflected and transmitted elec- 
tric field respectively for s and p polarizations. 



Fresnel coefficients for uniaxial material plates with 
a topological magnetoelectric response (axion) 

It is known from earlier work a^^i^^ that topological 
insulators contain, together with a dielectric response, 
a topological contribution to the magnetoelectric effect 
originated in an axion type Lagrangean. The inclusion of 
this coupling into the reflection coefficients was discussed 
by^ for the isotropic case. Their conclusion was that the 
tangential components of E and H were still conserved, 
although now H — fi~^B + aE where a is the topologi- 
cal magnetoelectric response or axion term and a is the 
fine structure constant (a — Analogously the inclu- 



sion of the new term to the boundary conditions given by 
(IC9)) . ((CT0)) .( (CTll and ([CT lt is i mplem ented by modify- 
ing the last two equations (jCll|) . (IC12l) (which come from 
the continuity of the tangential component of H) to give: 



(^11 -^ii): 



A, 



- R\\ 



{Ra 



A, 



hy + aCy, 



with hx given by (IC6p and hy given by 



(C14) 

(C15) 
(C16) 

(C17) 
It is now 



a matter of algebra to elucidate the refiection coefficients 
as in the last section. Inevitably, the relations become 
more messy. The reflection coefficients can be written in 
matrix form as: 



(r) 
s 

('■) 



2aq"fi±kz 



q"k,^i^c^ 2aq"fi^k, \(Ei'^\ 

{e^k, - q") + q') + q"kz,,^c? i E^f) ' ^^^^^ 



where Aa„ = (/x^fc^ + q') {e±k,_ + q") + q"kz^i±a'^, kl 



kl. q''- 



II 

and a = aO/n. 
To compute the Casimir energy, one should define these 



in the imaginary frequency axis and turn all the frequen- 
cies to a; = i^. Note as well that these reduce to the ordi- 
nary anisotropic coefficients presented in the last section 
when a ~ {). They also reduce to the isotropic reflection 
coefficients ([7]) when /x and e are isotropic. 
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